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A soicond order model for prediction of dissolved 
oxygen in a polluted stream has boon investigated in this 
study. This model uses Monod's ejcperession 3!or the consum- 
tion of organic matter, based on Micliaelir. Menton hypatlua- 
sis, to give the microbial concent rcatioii in a system, its 
due considerration. The resulting ecjuation has been solved 
numerically on IBM 7044/1401 digital computer system using 
Fortran IV language and the results compared vjith the exp- 
erimental data obtained on a simulated plug flow type str- 
eam. The values have also been calculated by the StreetGr- 
Phelps first order model for dissolved oxygen for comiaar- 
ison . . 

Theiil^dy shows that the model proposed: here ogives 
a better approximation to experimental data and hence sho-: 
uld better be used for accurate v/oric. 



Qj;' l;i jLi tJi-iOLjo 

t = ‘X'inie, lioar;j 

a = Subatrate coticencration at any iiistant, mg/i 

X = uiicrobial cioucentratiou at any inatant, uig/i 

Sq - Initial nubibtraub concentrauioii (at t=0) / n\g/i 

^ Initial iiiicroDial coucanci'ation (at t=^o) , iXig/l 

p. ==s GrowtVi rate coastant^ hr 

, -I 

p.^, i'ia:ivimuiu value ui the urowth race coiistaiit/ hr 
Kj^ = . i'iGchaoiis - hti'utoii constant 

K., = Substrate consurued per utiit iaasa of microblai mans 

synthesised or reciprocal of the yield coei'licieut/ ! 

= B.O.D* reaction rate, constant in Streeter - Phe Ip' s raono- 

— 1 ^ 

molecular fonmlation# hour 
K 2 = Reaeratiori rate cons taut/ hr”^ 

. I 

Ke = Bvidogenouo death rate, hr I 

K' = Ultimate oxygen demand per unit mass of substrate (glucose) 

metabolised, mg 02 /mg glucose as C.O»b* 

K" « Ultimate oxygen demand per unit mass of microbial mass 

metabolised eiidogewousiy, mg 02 /mgf microbial mass. 

G =s D.O. concentration at any instant, nig/l 4 

G = Saturation D.O. concentration under the existing physical p 

£> ■ i 

conditions, mg/1. 4|;; 

D = D»0, ciefiGit, mg/l. / | 
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II'iTRGOUi‘jTICl-1 

1.1 G;nt'JERM. ; 

The dicGolved oxygen cox'it'nnt. of a a-trearri, iP -a nutshell 
en its conditions and Uvsefulncssa froai the point of viev/ 
of human and animal consumption, industrial uses si.nd also from 
the viev/point of the ac^ntic life. Oxygen is a gns very npax"i- 
ngly soluble in v/ater, dissolving to the extent of only about 
8.4 mg/1 at nomal temperatures (25^0 at the saturation level. 
Upofi this meagre content of the die solved oxygen, depend the 
trillions and zillioiis of the forms of aerobic aipaatic lifve. 


present in a. stream, or anyotlier v;ater body, to keep themse- 
Ivos alive. Any substantial reduction, tberofore, in this vital 
ingradient to sustain aerobic aquatic life, already in short 
supply, can cause breaks in the natural food chain, resmlting 
in an upset in the othervjise existing balarxce of the mixed 
nature of the aquatic life. This may lead to e^rtinction by 
death or migrating of certain aquatic creatures ( 1)> which 
for ejrample, in the case of edible fish and crabs may mean 
an economic loss, in terms of their food value, for the comari- 
unity dwelling on the banks of the stream, A stream, lov? in 
its dissolved oxygen content also becomes unhealthy in its 
looks and givea rise to obnoxious odors. Its water becomes 
dark in colour due to the present of the products of . anaer- 
obic decompositioh of the organiG pollutants, and loses its 
Value, both recreational and as a source of water supply , 



The water can, oi' course, foe used for hunian consurription, after 
extensive treatment, but the cost factor becomes prohibitive 
and the corrirnunity is compelled to look for some other sources 
of supply. The property on the banks of such a stream/ may 
also lose its value, due to the odors from the streoiTi and its 
unhealthy look. 

In nature, clean waters are generally saturated v;ith 
the dissolved oxygen or nerirly so (2), clear and healthy to 
look at and v/ithout any odors* Vifhen domestic or industrial 
wastewaters are discharged into such water, a succession of 
changes in water quality takes place. The initial effect of 
the pollution is to degrade the physical quality of water, 
by the presence of dissolved or suspended pollutants. Gradu- 
ally, as the pollutants are utilized as food, and oxidised 
to gain energy for sustaining aquatic life, a shift to bio- 
chemical degradation takes place, and a lot of nev; organic 
compounds, some toxic to certain forms of aquatic life, are 
produced as a result of decomposition of the organic matter. 

In the biochemichal oxidation of the pollutants, the oxygen 
required is derived from the dissolved oxygen content of the 
stream. This causes a depletion in the D,0, content of the 
water. As soon as the D,0, falls below the saturation level, 
the atmospheric oxygen dissolves in water to make up for the 
loss. The surface layer becomes saturated almost instant- 
aneously and theirefore the rate of dissolution of the atmos- 
pheric oxygen in water, called reaeration, depends upon the 



tui.1:julenae prtsaent in it'/ wliich ietemiinn the pnvount of iT.iiting 
or renev/nl of tlie top layer in the stream, rioth the rate of 
deQ:r,ygenation of water uue to the biochemical degradation 
and the rate?, of reaerntion detorrrdne the amount of D,0, 
present in a stretch of polluted stream at any rnomont . Thus 
vie can picture .the? dis.solved oi-:ygcn .as an asset (or credit) 


and the biochemical o:-:ygen dematrd of the pollutants present 
as the liability (or debit) of the ntr-eam. By measuring tho 

D,0* content of a stream/ one can fairly judge, as to vdiat 

# 

eictent is the stream pjolluted and whether its v;ater can be 
economically put to the various beneficial uses. 

Thus./ the maintenance of a minimum 0,0, content 
provides us with a controlling criteri.a for the dis>charge 
of domestic and industrial wastewater in a stream. For 
optimisation of v/astev/atvcr treatment facilities and the 
optimum utilisation of the c-apacity of the stream to as,simi- 
late pol.lution/ v/i thou t any harmful effects/ it is nece.s.sary 
to predict v/ith a fair accuracy/ the D*0« concentratipn down- 
stream of a wastev;ater outfall. This ‘prediction has for long 
been based principally on the cl.assical first order bioche- 
mical deg radaition reaction kinetics theory forv;arded by Stree- 
ter and Phelps (3,4) , The Streeter Phelps oguation is based 
on the assumption that there are only two major processes 
taking place, which are: (a) consumption of oxygen in the 
satis faction of B ,0 .p . along the st retch by the biochemical p 
oxidation of the organic matter; and (b) repiacement of oxygen 
by reaeratibn dt the surface. Recently it has beerii pointed 
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out by many research v;oi3cers^ that the processes other than 
the tv;o assumed to be pc ©sent by Streeter and Phelps are 
also taking place sirriultaneously and must be accounted for 
while assesing the quality of v;ater in a stream. Some bf 
these processes are . (5) i 

(i) r6moval by sediniGinta’fcioa and/or adsorption. 

(ii) B.o.De addition by bottom deposit scour, and diff- 
usion of partly decomposed products of benthal 
decomposition . 

(iil) Addition of B.O.D. along the stretch by local 
runoff. 

(iv) Removal of oxygen from water to aerobic zone 
of benthal decomposition. 

(v) Removal: of oxygen from water by purging action 
of gasses rising from the benthal layer. 

(vl) Oxygen addition by photosynthesis in plankton. 

. '! • ■ ■ ■ 

(vii) Oxygen removal by respiration of plankton and 

fixed plants. 

(viii) Continuous redistribution of both the B.O.D. and 
oxygen by the effect of longitudinal mixing# 

The generally accepted model for predicting dissolved 
oxygen concentration in a polluted stream, given by equation 
2.1.3, developed by HW Streeter and EB Phelps (3,6,7), is based 
on the first order reactibn kinetics for both, the exertion of 
biochemical demand (8,9 ,10) and for reaeratlon < 6) ^ The use of 
a first order reaction kinetics for the exertion of B .O.D .has 
been questioned by a nu)rter of investigators, and alternative 



higher order reaction kinetics have been proposed (11,12,13,14), 
The process of toiochemical oxidation ot organic matter involves! 

(a) conversion of the substrate into biomass with varying degree 
of synthesis and (b) utilization of stored decomposition products 
and of cell structure* The later phase involves both the oxidative 
micro-organisnn and their predators. This clearly reflects on 
the relative importance of the concentration of micro-organisms 
too. Its importance can be better estimated by imagining the 
substrate available as the raw material in a factory and the 
micro-organisms as the workers, to convert the raw material 
into the finished products, carbon-di -oxide and water in the 
case of the stream. 

Let us imagine a sterile medium containing some bio- 
degradable matter, and enough D,0, to sustain aerobic life, had 
it been there* Since there are no micro-organisms or biochemical 
agencies present, to oxidise the substrate, no B.O.D, will be 
exerted and there would be no change in its D*o, content. Now, 
suppose, in two identical and equal portions of this medium, 
we introduce micro-organisms capable of degrading this subst- 
rate, one organism in one portion and two in the other. Natu- 
rally the substrate should be expected to be consumed and in 
the process, B,o,D, exerted, in the second sample at a rate, 
twice of that in the first one. Projecting the same to the 
normal concentrations of micro-organisms, present in natural 
waters, the importance of consideration of the concentration^^^^^^^^^^^^^^^^ 
of micro-organisras dan be appreciated. It becomes even more 
obvious when on© bears in mind th^t one of the factors aff ectihg 
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The D ,0 * in a v/atei- body is the exertion of B,04D.; avid fob 
wliicl'i vje should have a substrate to be degraded and the mici'- 


oorganisms to do it, 

A second order raodel, involving both, the ooncentrati- 
ons of substiaite and that of the micro-organisms, is proijosed 


and investigated in this study. 


t takes into consideration the 


microbial concentration, according bo theMonod's expression 
for consumption of organic matter, based upon Mechaelis-Henton 
I’lypotliesis » Monod’s work (15) on microbial giroeths in Chcmosbats 
has contributed significantly tov.-ards a deeper understanding of 
tViO principle.s involved, and the role of the concentration of the 
micro-organisras in the exertion of B.O.D, Recent experiments on 
tVie study of concentration effects in the biological oiciclation 
of trade v/astes by VJilson (16) have extended honod's results 
obtained for purs cultures, to heterogenous mixed cultures, 
thus making them apjjilicablc to Sanitary Sng inhering problems. 

The term representing the dissolved oxygen consumed in 
the biochemicral oxidation of the substrate is derived from the 


Monod’s expression, v?hicVi involves the concentration of both; 
the substrate and the micro-organisms. Thus the effect of the 
microbial population on the dissolved oxygen is accounted for. 
As the biochemical degradation of even sinipic compounds like 
glucose consists of many intermediate steps, involving numerous 
intenaediate products, with their ovjn rate of reaction the 
slov^est of which becomer. the governing criteria for the overall 
rate of reaction » But t o dote rmine the exact rate o £ reaction , 
one will first have to idehtify this: slowest step and measure 
the concentrations of varipus enzymes involved, react ants an 
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intermediate products at that stage. Kov;/ this may not always 
be an easy job, specially in the light of the fact that the 
whole chain of reactions, the. slowest step and the involved 
chemical and biochemical agencies will be different for the 
different pairs of substrate and rnicro-ox'ganisms. These 
complications will further be increased in the case of mixed 
cultures of micro-organisms, involving various types and spe- 
cies of organisms. The situation will become even more comp- 
lex for the complex substrates like proteins and other nat- 
ural products. If we make a reasonable assumption, that the 
microbial concentration represents the overall biochemical 
activity, which is justified, because afterall, all the enzy- 
mes and, other necessary biochemical agencies have to come 
from the micro-organisms themselves, then we can safely say 

that the rate of reaction will depend upon the concentration 

♦ 

of substrate and the micro-organisms present at any instant 
and hence, the reaction will be a second order one. 

1.2 AIM AND SCOPS OF THE STUDY 

The main aim of the author in taking up this study 
has been to investigate into the effect of both the microbial 
population and the substrate concentration on the dissolved 
oxygen content of a polluted stream* A comparison of this 
rational model with the Streeter - Phelps equation is also 
to be made, with the primary aim of bringing out the inadeq- 
uacy of the later formulation in dealing with the practical 
situations..; ^ 

This study has been limited maihly to the following 



horizons: 


1. To formulate a rational model for dissolved 
oxygen profile in a stream/ giving /the microbial 
concentration its due recognition, 

2 , To verify the model experimentally on a simulated 
model of a plug flov; type stream. 

The numerical solutions of the mathematical formula- 
tions are worlced out on the IBM 7044 digital computer , A 
fev// model dissolved oxygen profiles for varying initial 
conditions are plotted experimentally and compared with 
the theoretical solutions of the proposed model for verifi- 
cation. 



CHAPTER II 


LITERATURE REVIEW 


2 , 1 FIRST ORDBEi KIRETICS FOR B.O.D, REACTION £c BQUATIOM 
FOR D.O, PROFILE BASED OH THAT: 


The classical/ first order theory for the reaction 
kinetics of the exertion of B.o.D,, developed by Streeter 
and Phelps (3,4) states: "The rate of the biochemical oxid- 
ation of organic matter is profpiortional to the remaining con- 
centration of unoxidised substances# measured in terms of oxi- 
dizability" (4) . This law, expressed in the differential form 


becomes : 


cX L 


or - 


where, l a B,0»D* present at any instant t . 

Constant of proportionality, the reaction 

rate constant 

Integrating from, t » O, L = L_. to t » t, L = L, we get 


L — ® 


-kit 




The expression representing the dissolved oxygen profile 
(or sag curve) , based on this theory, developed by H,W« Streeter 
(6,3) and, analysed by G.M, Pair (7) is given by 






.2,1.2 


which upon integration yields , 


D = + . . 2.1,3 
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V7here, D « Oxygen deficit, at a point distant, time of 
flow t, 

Dg~ Oxygen deficit at the reference point, at t = o 
Ultimate first stage B.O*D« at the reference 
point, at t - o 

}C2® reaeration rate constant 


2 « 2 EFFSOr OF SOOURIMO, PJSPOSXTIOM, PHOTOSYNTHESIS, AND 
BENTHAL DECOMPOSITION ; 

H*A. Thomas, Jr», (21) in 1948, pointed out that, 
deposition and scouring of the suspended matter and stream 
bed sediment respectively, must also be taken into account, 
and included another constant k^, in equation 2,1,3 above. He 
postulated that the rate of deposition and scouring is 
proportional to the remaining B.O.D. and modified equation 2,1,3 
to the form: 


k2~(l^l+k3) 





2 , 2,1 


where, k 3 « the constant of proportionality, reflecting the 
composition of waste and the receiving Water 
(regarding the suspended matter) , and the quiesence 
of the streanf at the point under consideration. 

The negative value of k^ , indicated scour, reversing the 
process "ot depositlon,, v:V 

( 5) brought up the importance of addition 
of B,0 .D * from benthal la]^r and local runoff f ^ 

Oxygen by benthal demand and ; plant respliration ; and addition r 
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of oxygen by photosynthesis. He made the follov/ing assumptions: 

1. The Stream flow is steady and uniform. 

2. The process for the stretch as a vjhole is a steady 
state process/ the conditions at every cross sec- 
tion being unchanged V7ith time. 

3. The removal of D.O.D. by both the bacterial oxid- 
ation and the sedimentation or adsorption or both, 
are first order reactions/ the rates of removal 

at any section being proportional to the amount 
present, 

4. The remoVfal of oxygen by benthal demand and by 
plant respiration, the addition of oxygen by pho- 
tosynthesis and the addition of D.O.D, from the 
benthal layer or the local runoff are all uniform 
along the stretch. 

5. The B.O.D, and oxygen are uniformly distributed 
over each cross section/ thus permitting the 
equation to be written in the usual one dimen-* 
sional form. 

Based on the above assumptions Dobbins proposed the 
following expression for the prediction of B.O.D* remaining 
and oxygen deficit down stream a pollution point: 




k2"' (Ic^+Kj 


HsCKl+kg) 


) ( ^ ) # • % ♦ 2 # 2 * 2 • 1 
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L«L^e ^ ) ••« 2. 2. 2. 2 


v/here, L « B.O,D. remaining at a time t 
hj^~ Initial B.O.D, 

rate of addition of B.O.D. along the stretch 
D = Dissolved oxygen saturation Deficit 
Dq~ Initial Dissolved Oxygen Saturation Deficit 
Pj 3 “ Rate of removal of D«0» along the stretch by 
benthal demand and the effect of plants 
k^= F<ate of constant for trological oxidation 
k. 2 ~ Reaeration constant 

kj* rate constant for B.O.D. removal bysedimentation 
or adsorption. 

Another equation for the D.o. profile between the 
stations a and b proposed by Camp (22) in 1965, includes 
and accentuates the photo synthetic reaeration of water 
due to algae and other aquatic plants/ as compared to the 
atmospheric reaeration/ is 


” k2~()ci4k3!r ^ 2.3(ki+k3) ? ( 10 ^ 


t ( 


- 7r4“r ^ ^ , . 2 .2 .3 

^l+kaT 2.3 ki ^ 
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" Average oxygen saturation clGficit# over the 

cross section at river station 'a', Steady state, 
in ppm 

= Average oxygen saturation deficit, over the cross 
section at river station steady state, in ppm 

= deoxygenation constant, day"”^ 
k 2 “ atmospheric reaeration constant, day*"^ 

rate constant for settling out of B.O,D.^day 
I'a =* average steady state first stage oxygen demand 
at river station a, ppm 

p = rate of addition of B,o»D# to overlying water 
from bottom deposites, ppm/day 
c< a rate of production of D.o, by algae through 
photosynthesis, 

U e average velocity of river between stations a 
and b, miles/day 

^ « longitudinal mixing coefficient between river 
stations ‘a* and 'b', square miles/day 
X = distance between stations 'a' and Vb', miles 
0‘Uonnel et*al. (23), proposed the following model 
for the D«0« profile, based upon the studies of the Truckee 
River in Nevada State, U,S.A», stressing upon the effect of 
benthic algae oh the stream D,o« 


k j_lt 






2,2.4 




1 4 


vi7here/ P~R = net ojcygen contribution toy i:)hoto synthetic 
activity 

P = oxygen gain toy x>hotOGynth03is 
R = oxygen consuniGd in algal respiration 

dissolved oxygen saturation level deficit 
D,^= Initial dissolved oxygan saturation level 
deficit 

kj,= deoxygenation constant 
and/ k 2 = reaeration constant 


Akerlindh (24) emphasised the importance of the 
importance of the immediate oxygen demand and the oxygen 
demand exerted toy the bottom sludge/ and proposed the 
following equation for the D.O, profile 


‘r 




i-'J (e 


-kit 


” ® ) + k 


D. . . . . 2.2. 

kr-1% ® 


(• 

!• 


v/here, » Initial dissolved oxygen saturation deficit/ ppm 

Lj^ = Initial immediate oxygen demand, ppm 

= Initial B«O.D, in the stream/ due to “flowing 
pollutional load'*/ i*e. dissolved matter and 
non-r.deposi ting suspended matter, ppm ^ 

= B,O.D. in the stream due to ''depositing pollut- 
^^;ional':load*'/■ppm:,^^■■V■;■■'.■ 
kj|_ != specific fate coef f idient of dsoxygenatioh due to L 
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= specific rate coefficient of deoxygenation 

due to L, 
f> 

kg = specific rate coefficient of deoxygenation 

due to L 

s 

k^ =1 specific rate coefficient of reaex'ation 
t = time of flov7 in days, 

2,3 SECOMD OR HIGHER OHDER E;g>RSS;310MS FOR B,O.D. 

RBACJTIOM KIMETXCS : 

Many research - workers have questioned Streeter - 
Phelps first order B,O.D, reaction theory/ in order to arrive 
at a tetter coordination of theoretical and experimental data 
of exertion of B.O,D, and prediction of D,0, profile. Several 
alternative expressions have been proposed/ some of which are 
alternative expression for the exesrtion of B,o,D. as 

,y^ *S(Mlogt'fB) • • fl* ,0 ♦ • o 2,3,1 

v/here, * ultimate B ,0,0# at the end of the stage 

S =s 5~day B.O.D, intercept of the line 

.. m ■ i 

M « = B,0«D, rate parameter 

Vs ■ ' 

B « — = B,0.D. rate parameter 1 

S’ 

m ss slope of the plot of v/s log t (straight line) | 

b » intercept of the plot of v/s log t t 

The values of M and B were found to be 0,85 and 0,41 respect- [ 
ively for the domestic sewage* /I 

Suggesting that the rate of removal of sustrate is ^ | 
propbrtionai to the Square of the remaining substrate/ Young (13)| 
proposed the following second order differeriitial equation fdr | 
the axertibn^pf . ::r.| 
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^ ^ “ Y ) - ^ ~ y) , ^ ^ ^ 2,3.2 

Integrating, w© get 



where, L a initial substrate concentration Cat t =* o) 
t a time 

y a B.O.D. satisfied by time t 

ReVelle et.al. (14) emphasised on the importance of; 
(i) The rate of removal of the substrate is proporti- 
onal to the product of the concentration of subs- 
trate and bacteria, 

(il) The bacterial population is proportional to the 
amount of food removed. 

On the basis of these assumptions# they proposed the follow- 
ing e>^ression for the exertion of B,O.D, 



b + l> 

b 3k(L+b)t^j 


, 2.3.4 


where, y <= B.O.D. €:<p res sed at time t (days); mg/1. 

L = Ultimate 3.0. b.# mg/1 

b = a constant representative of initial bacterial 
populetipn# mg/l 
k = rate constant ( (mg/l) day) 
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2.4 SECOKD ORDKR D.O. PROFILE EQUi'.TXQlNS : 

Dresnack et.al. (Jfl) the £ollov;ing tv/o 

second order partial differential equations for B.O.D. 
and D.O. along a stretch of polluted natural stream/ 
analogue to the equation representing the conduction 
of heat in solids: 


ai 

at' 

= D, . 

« uah . 
ax 

- D+Lq , . . 2. 4, 1.1 

a_c 

. Dj, . 

“ - 

' k,L -I- k^ (C^~C)^Dp- , 2. 4. 1.2 

dt 

" X 

fax 



where, L as first stage B.O.D., ppm 

C =s dissolved oxygen concentration, ppm 

Saturation dissolved oxygen concentration, ppm 
coefficient of longitudinal dispersion, aqri-iiles/day 
Uk average stream velocity, raxles/day 

-1 

kj^sa B.O.D. reaction rate constant, day 
k^- reaeration rate constant, day 

k “ rate of removal of B.O.D. by sedimentation and 

3 ■ ' 

-1 ■ 

adsorption or both, day 

L = rate of addition of B.O.D. along the stretch 
a 

D] 3 = Net rate of removal of D.O, by all processes 

other than biochemical oxidation of the flowing 

X = distance along the stretch, miles . 

/ /^t; «;:tirne/ /:da:^ 

Equations 14 and 15, aolyed difference method, ; i 

yield:;'" ■ " 





Li^ , M''! ( 


k(ki+>c3) 


• k D. 




kdcj^^k^) 2]cDi, 




■I* ^i-l/n + k L^ . . , 2, 4. 1.3 


and 


(1+ -^2) „ Cj^-1, n n (1~ 


k k2 2lrJDr 
— 5^1 

ir 


-KJj_+l, n(|^^)+k k2 Cg-k DB-(ti/n'j-Lj_,ti+l) , ,2. 4, 1.4 
v;hex'e# k = time mesh interval 

2,5 OTHER EXPRESSIONS FOR D.O. PROFILE ; 

Many workers have avoided the usual deoxygenation ~ 

* ' ^ 

reoxygenation approach in giving suitable expressions for 

the dissolved oxygen profiles and gave some emperical 

■ ■ ■ . 

formulae or some expressions based on some other" kinds 

■ ^ I 

of theory. An expression proposed by Churchill and Buck- 
ingham (25)# thus completely avoids the use of the usual 
deoxygenation and reoxygenation coefficients , k^ and kj^# 
by making a statisfical fit of observed data according to 
the assumed relations: 

D,0. drop * a + b L(5) tCT + d/^ . , * , 2.5,1 [ 

where D ,0 , drop = difference in ppm# between the P ,0 . at I 

a Igcatipn above the single pollution i 

source and the D.p, level at the criticr' i 
':- . :al, point' of :;the sag; ' cu'rve ■ V 
L(5) - Five day B.Q.p, in ppm# measured at a . ^ 

• location in the vicinity of the , critical , 
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T s= Water temperature, ocj 
Q = Streani flovj, cEs 

a,b,c,d != constants to be determined by statistical 

analysis of data obtained from a set of obser"ved 
D,C, profiles. 

For a stream v/ith a discharge that is steady at 
any particular cross section/ but v;hich may vary along 
the course due to the distributed surface runoff and 
groundv/ater# Wen-Hsiung hi ( 26 ) gave the following expression: 

C3-C = e-’^2[c3-J'C5 ) + | 1 CO^-Cj,) 

& 

X X d 

+fl5> ) 7)^1 +J ^ ’‘1 ^ ® . . . 2.5.2 


Where/ =D.O. at saturation 

C(x,t)= D.O, in mass per unit volume of water 
^a(x) 2= D.O, of added discharge 

F(t) = D.O, of stream water at out fall i.e. C(o,t) 
f(t) = B.O.D, at cross section of out fall after 


miximg with stream. 

Y ? dx 

) - t > J » t-T 

t “time 

X 5= distance along the stream 
V{x) =» speed of flow = Q/A 
A(x) « Cross sectional area of stream 

K^Cx)* Dilution coefficient ^ 

. :-‘A: dx 
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V 


dx 




V 


dx 


^ -2^ dx - ! 

/ 

This model assums that the D,o, and B.O.D. of added/ 

f 

surcharge along the stream, downsti'eam to the outfall con 9 'e~ 
rned is steady at each point, but not necessarily unifor^. 
along the course of the stream* The D,o, and B.O.D. at the 
outfall may fluctuate with time. Well mixed conditions, both 
vertically and laterally have been assumed, so that the 
problera is essentially one dimensional , 

Thoman (27) has proposed the following mathematical 
model based on the concepts of “systems analysis" for the 
D,0. in streams 

C(t) = C^g A(w) Gos(Wt + © M(w) ) . . . 2. 5*3.1 

This is a linear system, much oversimplified and 
is rarely met in nature ♦ Thoman goes ahead to generalise 
this model to the form 

c (t) =a£jok[&xi'«'loLhij(w)]i,j,i;os[wt+|Bki(w)}ji, 

J. ctj . : 

+ ^©ij(w) j . . . . 2.5*3,2 

where, Gj(t) = p*o. in m 

“ Saturation D*G. varying in a periodic form 
with the temperature of water, over the year 

‘ « angular fteguehey of 

©^j(wX « phase ahgie shift of Cg^{t> 
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Aj^jCv/) a Arftplitude ^attenuation oJ: fj. (t) for 
solution of Cj(t) 

d^(t> = Decay coefficient in system i 

^^(t) =5 treatment plant forcing function in system i 


A simple emperical equation of calculation of D,0 
drop dov/nstream to a pollution, vms given by Smith et.al 
as follows: 


Y 


K Q 2 


T/10 


2*5,4 


where, Y = D,o, drop in ppm 

B ss B.O.D. loading in thousands of Ibs/day 
F =s Stream flov; in mg/d 
T as Temperature of stream in 
K = proportionality factor 


(2B) 


Okun et.al (29) studied the waters of Catawba 
river in south Carolina, O .s«A« and fitted the 188 claya 
recorded, to a general equation given below, and obtained 
a multiple correlation coefficient of 0,970. 

Y » 9.47-0.49 07 Xj_X]^+0,0i9 X 1 X 3 

-7,1 Xj^Xjj-0. 341 x^ -0.00035 X 2 X 2 +0. 00528 X2X^4, 3 x 2 X 4 
-0.0839 X 3 +0.000 17 X 3 X 3 - 2 .O X 3 X 3 - 2.0 X^X^+SO X 4 

^ ^ ^ . . , • ; 2 . 5 . 

Where f Y »■ D.0.'’ln\mg/i; 

Xj^ss. vIntake^;P:.6 .vlmg/l ; ' : ':v ■ j'v; ■/ 

■■ 352 =®' "River ■ flw I 000 ::cfs ^ 

X 3 « Temperature/:, ^ /-.v'//’': S- 

'■X 4 S!; /Discharge; ;r At. ip - West e. 



n *7 

tr j £-• 


This equation "being mathematically cumbersome/ v/as farther 
simplified by the authors/ eliminating several terms which 
have little value for prediction purposes without any loss 
•of accuracy, and its final form came to be 

Y =s 8.88 - 0.44 + 0.068 x^XjtO.Ol2 Xj^X2~0.19 X 2 

- 0.059 X3~32.7 x^ o . * ,2, 5. 5. 2 

O'Connor (30) v;ith the help of a set of linear 
differential equations, evolved the follov/ing equations 
for prediction of oxygen deficit in the streams: 

^ £ ^e”^'^Sinh(jd)] - ^|2®~^^^Sinh • 2.5.6. 

for - 00 ^ X C -d and +d 5 . x < + ^ 

and 

D = ^ (Jx)] - 

for - d ^ X £s + d 
where, D = oxygen deficit 

. : ■ ■ ■ '+jX 

1 = B.O.D. concentration - la ® 

W .' 

.la“: 2 A.jKE. ^ 

■ , , W ■ 

■ ^6.^ ?sr ;■■■ 

v' ; '-j == K/E::= fluX' , 

W = daily mass rate of discharge of B.O.D, 

A : « GroiSs ^sectional area:- ; : , , v-; : 

E » ::dispersie>n : coefficient 



” reaeration coefficient 
X = distance along the stream 
Uniform load is assumed to be from ~ d to + d* 
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CHAPTER III 
THEORETICAL BACKGROUND 


The model investigated in this study is based on 
Konod’s expression describing the kinetics of beooxidation 
of organic matter as described belov/# 


Defining ^ as the rate of grov/th of nucroorganisms 
in lograthmic phase/ or 

i*= is 

v;here/ x = concent x'ation of microorganisms 
^ = grovjth rate constant 

we have/ 

1 dx 

P ^ dt ^ 

dX , 

OX/ X • • * # ► • * • 3 »1 

dt “ 


Monod (15) shovjed/ that, for pure culturas, which 


was later extended by vJilson ( 16) to mixed cultures, appli- 
cable to sanitary engineering systems. 



Hmax ^ 

K„ + 


3*2 


where, u » maximum value of grov/th rate constant (W) 

v/ith uni imitted substrate* 

Kg = Michaelis - Menton const ant, given by subs t~ 
rate concehtrabion at which the observed 
growth rata is ' one half of its maximum 



Goinbining the equations 3.1 and 3.2/ v/e get 


ax 


„ Mmax S 

U'«. 'r o 




. . . 3.3 


Also since. 


X 


^ K:: 


s) 


3.4 


v/here, ~ initial concentration of rfsicroorganisrns 

* initial concentration of substx'ate 
K^, is substrate consuraed per unit mass of micro- 
organisms synthesised, or teciprocai of 
coefficient y , ", 

X/ i3 

Expressing in a diffox'ential form equation 3.4 becomes 


^ „ v dx 

dt dt 


3.5 


Substituting the value of dx/dt from equation 3.3, and 
the value of x from equation 3,4, we get 

S{A-S) 


ds 

dt 




max q .«3 


♦ ♦ ♦ • 


. . . 3.6 


where. 


A = Xo + So 


Integrating* equation 3 ,6, we get 


t 


V f , A-S ^ A.,s 1 


;■ . /x :.';, , 3 ..7 


*S t ep s o f int ration sho^ in the app en<ai x ; 
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YHB DISSOLVED OXYG^W MCDi;::L. : 
vvriting a rnas.-j 'balanc<i oguation, for the dissolved 
oxygen, for an aerobic biocheinicai nyGtC;in, v;e have. 


Cliange in 
oxygen content 
of the cystem . 


Ci-iange in the oxygen I Change in the oxygen 
content due to subs- -i- content due to end- 
trate metabolism J Logeour. metabolism 

Change in the oxygen 
t content due to x'eae-- 
ration from atmosphors 


^dt^ovorall” ^ ^dt^ substrate ''' ^ endogenous 

metabolism metabolism 


^ reaerat 


atxon 


where, K, 


dC dS dC dx t V 

K2 (Cg - C) 

dS dt dx dt 


reacration rate constant, day 

dissolved oxygen concentration at time t , mg/l , 
saturation dissolved oxygen concentration 
under the existing physical conditions, mg/l. 


As under endogenous conditions (31) , 




. V ... '. 3.1.2, 


where, K 


Endogenous death rate constant 


Substituting the values of ^ from equation 3 « 6 and that of 
from equation 3 . 1.2 ip equation 3 • 1,1, we get , 
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dt 


SCA-S) dC 
rrnax p — TS — • 


''® §; ■‘' ‘-2 *''3-^^ 


• • * ♦ • 3*1*3 

t'i'ow/ substituting -. C - D, the dissolved ox'ygen de-ficit 


vje getf 



- 

dD 

dt 

=5 „ 

j^rna;j 

S(A~S) 

/ ~ K . X 

dc 

dt ^ , st*^ 

■ _ "1" ^‘••ry D 

d3C “ 


dt 

+ Ko 



S(A-S) 

dC , 

dC 

or, 

D 

Pm ax 

K«+ S 

ds 

dx 


dO 


D = 


C! V. 

u> \ rv.*r?iM>y 


i (s,j-s)j ^ 

Kx y dx 

or, 

dt 

+ K2 

Pm ax 

Iviri 

K' + Ke 

or, 

dD 

dt 

+ K2 

D == 

pm ax 

K* :i?(A-.S) 

. + he^ (A^S) 







Kx 





3 . 1*4 


v;here, K' 


as 

ac 


- ultimate oxygen demand per unit mass 
of substrate. 


K*‘ = H;- = ultimate oxygen demand pel* unit mass 

dX 

of microbirganisrfts under endogenous 
coxiditions. 

The above equation can also be written as 


or^ 


fiS , ^ + %,2 D 

dS : dS/et • 


^ pmax ■ 


-s : 

icT^ 


Kg K” 

: Kx:-: 


"7 

J 


- -n ^ W :k^ 1 v.;a..S ' • 



2a 


sul^oti'cuting tliG vtiluG o'f riTom eqjuati-on 3»5 




Pmax 


k_+ 3 ^ i ; 


,nK” ? 

rcj-. \ 


'Pmax 

*‘>‘^4' >«■) 


® , B . ?K' . .] 

[*raax SIA-SJ C K:: J 




3.1.5 


This equation has been solved numerically/ as the analytical 
solution becomes quite complicated and curabersimvj for the 
purpose of obtaining the data . It v/as solvced numerically 
by Rungc-kutta fourth order method on IBM 7044/1401 digital 
computer system using Portain IV language. 


DETBt^^lRATlON OF THE a0iMf3T/VTT S 

Equation 3.3 upon integration yields (30), 


In X - prnax t - In \ 


where, A = Xq + Sq 


3, 2. 1,1 


Rearranging this equation, v?e get. 


1 , /X V 


Pmax 


Ka/A i 


1+Ks/A l+Kg/A t 


• In 


. • »; .' 3»2»1»;2. 


Now substituting, . 

■■::"^;::'.v:;;:;:b, k^so 


h «i x x^ : ■ 

pmax/^l+^sA) 


, J. . .V v: :3;;2 »l .2 V a ; 
/■ ■, ' c .V ■ ' v; 1 ‘ d . 2 .1 *2: b : 

♦' J.; • * ' V 3 ' . 2 . 1 ;* 2 ' .-C 



n 


l+KsA 


3. 2*1.2 d 


equation 3, 2. 1.2 becomes 

In h:/z.:o) = n ^ in (l~bh) +m 


3. 2. 1,2 e 


Tills equation is 
and v;ill plot as 


linear in ' In C:v'>to) and i In (1-bh) 

U 

a straight line for there two variable 


grouiJS* The value of b can not be determined by direct mea'~ 
surement/ and has to be fixed by a trial and error procedure. 
For this In (x/xh^) and i in (l-bh) are plotted on an 


arithmatic graph paper (Fig. 1 ) for various values of b, 
till a straight line with a positive sloi^e is obtained. 

Then using this value of b, the various constants are 
calculated from the follo'wing ecjuations* 

m 

IHnaK • • • • 3. 2.1.2 f 

K _ s: * * « * • • » • • ® * 3. 2. 1.2 g 

s l-n 

■ Iv b 3.2il.2 li ' 

3.2.2 K* . and K’* 

Since the second terra of the equation 3 .1,1, repre- 
senting the oxygen consumed in endogeous metabolism is^^ ^ 

{§ 0 .} = dx 

dt endogeous "" av * fH- ^ 

, metabolism ; 

. — : Kq K.**- X.: , • . . ; . , /» ■: * 3 . 2 * 2 » 1 : . 

^A1 so because, dC:.i= iP.-r. v ; r;-;-?, 

, ■ -dxpP'r'K^b'bp 

we can conclude that '.:k" " x 



sufeGtituting in equation 3. 2. 2.1 vje get. 


Therefore, 

dc 

dt 


i\Q C 


3.2 .2 


n 


v/liicVi is a first order kinetics equation v/ith reaction 
rate constant as . Tlierefore the constant can be 
determined by the Fa jinioto' s Method (17) . 

Since K*' is defined as 


^ " ultiimate oxygen demand of the microbial 

mass under endogenous conditions. 

Therefore K” v.’ill foe given by the ultimate oxygen demand 
ho, in the Fu jirnoto* s Method. 


3.2.3 B.O.D. RSAoTIOK RATS COhSTTM^iT 

According to equation 2.1.1, the first order. expre- 
ssion for the kinetics for the B.O.D. reaction is 


L == L- e 

rJL 


-kit 


Taking logi'atbins 

In ” *” ^1 <1 ^ ... « . 3 .2 .3 .l a 


As L represents the B.O.D. remaining at any time t, which 
is representative of the remaining substrate concentration. 
Therefore replacing L by S and by be consistent with 
the terminology followed so far, we get 


in (S/Sq) « - t 


3.2 .3.1 b 


Thus a straight line will be obtained by plotting In (S/Sq) 
against t whoso slop® will give 
rate const ant 4 
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3,2,4 i;SASRATXOK C OH ?^T /^r4T, Kp 


The rate of reaei'ation or x'ate of change 


in dissolved 


ojcygen concentration or oxygen deficit, is proportional to 
the oxygen deficit cjvisting at any time t (2) , Expressed 
mathematically, 

~ Ko D 3, 2 ,4,1 

dt 


Ihtegrating between limiLs from; at t =0, D - to at 
t - t, D = D v/e get 

In {,D/l)q/ ® ■" ^2 t *«>,,* 3,2 ,4*2 a. 

.or D = e-*^2 t 3.2.4,2 fo , 

v;here, D « - c =; D,0. deficit at any time t, rng/1 

^'S "* Cq= D.O, deficit at time t = O, rag/l 
Cg= Saturation D,0, level under the given physical 
conditions, mg/1 

c w D,0, concentration at time t, mg/1 

Cq- initial D,o. concentration at time t =0, mg/1 

Kgs constant of proportionedity, the reaeration 

^ 1 ' 

rate constant, day f 

The data of dissolved oxygen at regular intervals, in 
deoxygenated water, when i^lotted as In ~ agaist t, will yield 

■WQ . 

a straight line, vAiose slope, will give ^ ^ 



CWvPTER IV 


MATERIALS AI<ID METHODS 
4,1 EXPERIMENTAL TECI-miQUES At^D EQUIPMERT 








A solution of substrate Cgl\icoGe ) , c.O.D. 400 mg/1 
v;as taken v;ith k 2 HP 0 ^-idi^F 0 ^ buffer and ammonium sulfate 
solution to provide phosphorus# control of pll and nitrogen 
for synthesis# v;as taken in 30 cm*r 10 era. 0 cylindrical jars, 
and inoculated v;ith knwvn araount of seed microorganisms# pi'c-. ■ 
pared as described above in the section ^•.1.2, This v/a.s sti- 
rred by a flocculation paddle stirrer# to keep the v/hole 
thing mixed and aerated. The concentration of substrate as 
c;,0,D, and that of the microorgvanisms v;as determined at 
regular intea-vals of time, as decribed belovr in section 4,2.1 and 
4,2.2. This data vjas utilized to calculated the .constants 
f*max^ described in sections 3 . 2. 1 and 3,2*2 (Fig Jl ) , 

4. 1,3,1 K*^# K p. and K ” , 

The constants K' # k”. and Kg, were determined vjith the 
help of ivarburg Respirometer. The flask arrangement was kept 
as follows; 



Flask 

1 Flask 2 

(endogen- 

Flask 3 
(Thermoba 
ometer) 

1* Main compartment. 

ml 

ml 

ml 

microbial auspensioh 

1.5 

r:,;' ''::i.5 


distilled wat.er 

'':V- : 


2.5 

2 * Side ArTti/ glucose ^plutipn 



■■'■b 

3 , ' Central 'S^eil, 6M KOH solution 

0.1 

' Q.i,’'.'; 

0.1 




S , SBSTRAIS CONCEHXHASIOH aB C . 0 . B . , ffli 
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The respiration data obtained v/ith flask 1 was 
used to detentiine K' , the ultinicite ox^*gen uptake jjer unit 
mass of the substrate (glucose) / by I?ujimoto‘s Method. The 
product Kg and K'* v/as also determined by Pujimoto's Metliod, 
from tlie data obtained with flask 2 (Fig, I ) , 

, 4. 1.3, 3 RSAFRATXOK RATE CQKST.;\NT. Ko 
Uater v/as deoicygenated by adding sodium sulfite 
solution/ 7.88 mg/1 per mg/1 of dissolved oxygen to be 
removed/ under the catalytic of cobalt ion at a concert ra- 
tion 0.01 ing/1 (18) , Tills deoxygonated water taken in a 
2.5 R/ 30 cm, jc 10 cm. / cylindrical jar and aerated v.’ith the 
floculating stirrer at a known constant speed. The dissolved 
oxygen concentration was noted at regular inten/als of time 
and K 2 V^as calculated according to the equation 3,2.4,2a (Fig ,3) , 
4.1.4 THE SAG caRVS8 

A glucose solution of concentration of 25 mg/l vjas 
prepared in aerated v;ater/ and taken in five 30 cm x 10 cm / 
cylindrical jars/ to simulate a plug flow stream. These jars 
were inoculated with varying amounts of micro^jial suspension 
to obtain different S/X ratios. The initial microHal epneen- ■ 
t rat ion was s 

Jar 1 Jar 2 Jar 3 Jar 4 Jar 5 
Microbialconcentrati- 

2,1 4,2 6,3 8.4 10.5 

The D.O, concentration was measured as described 
section 4.2,3 for all these jars a the regular interval s 
of time.. The contents, of jars were aerated by aerated, by 
''a '.flocculating stirrer at 70 
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4,2 MTALYTXGAL TECI-RIICUBS : 

4.2.1 MXCROBI.AL COMCEMTRATXOM 

The gravimetric determination of microbial GOi'jcon- 
tration for stcindarctisation purposes v;as made by filtering 
a know volume of the microbial suspension under vacuum vrith 
the help of a filter pump, through a V/hatrnan - 42 filter paper. 
The residue on the already v/eighed filter paper was dried at 
120 ®C for two hours, cooled in a dessicator for 15 minutes 
and weighed again to determine the dry v^eight of biomass 
caught on the filter paper* The weight reading was taken 
after 5 minutes of taking the filter paper out of the dessic- 
ator, to standaridise the effect of moisture absorbed by the 
filter paper during the process of weighing. A calibration 
curve (Fig 4 ij ) of absorbance v/s microbial concentration as 
dry weight/lv determined gravimetrically as described above 
was prepared. The absorbance was meansured on a "spectronic - 
20" spectrophotometer*, at 400 mji wavelength. Later, during 
the e^jperiments, the mixed liquor of the jars was pipetted 
out and its absorbance measured on the spectrophotometer. 

The corresponding concentration was then read off the calib*- 
ration curve (Fig.^), 

4.2.2 SUBSTRATE COKCENTKATIQN AS G,Q.D . 

To determine the C.O.D, of the Jar liquor, to calcu- 
late the concentration of substrate remaining at any instant, 
the liquor was centri fuged at 5000 rpm for lO minutes , The 
of clear supernata was determined as per 
Methods of Examination of Water : and Waste Water" <^19 ) . 



^^Eausch,"& Lomb, U .S.A, - 
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4,2.3 DISSOLVaP MEAGUREMKHTS 

The dissolved oxygen concentration v/as measured by 
the help of a "Precision Galvanic Cell Oxygen Analyser"*, 

A .calibration curve (Fig, 5* ) of the instrument reading 
v/s the actual dissolved oxygen concentration as detemined 
by iodometric titx-ation method (20) was prepared. Later, 
during the experiments the instrument reading v/as noted 
from time to time and tlie corresponding dissolved oxygen 
concentration read off the calibration curve. 



^Precision Scientific : Co 


,, Chicago, U 








CHAPTER V 


RESULTS mP DISCUSSIONS 
5 . 1 BVALUATXOM OF THE C0MST7-:tjTS 

ij’rom the data ot the variation ot substrate and 

the iTiicroDial concentration, obtained as desci'ibed in section 

4, 2.. 3.1, In (S/Sq) versus t is plotted in Figure Ho. ^ as 

a straight time. The B.o.D. reaction rate constvant as 

described in section 3,2,3 is given by the slope oi this 

line to be 0,04606 hour*"^. For calculating w„.„, and K„, 

1 1 

± In (x/Xq) and r (l-bh) are plotted on y and x axes 
fc ^ t 

'respectively as described in sefction 3,2.1, for various values 
of b, till a straight line vjith a positive slope is obtained 
(Figure "? ) . This value of b comes out to bo p.OOG, Hov/ 

Kj 3 and Kjj. are calculated by the- equations 3 ,2 .1,2 f, 3. 2,1,2 g 
and 3. 2. 1,3 h and the following values are abtained, 

Pmax ~ hr*"^ ~ 15.6 day"'^ 

tt'; ; Kg ■■^ ■415, o'mg/r .'' 

'; 2 ,64; 

' ■ . : 5,1.2 K^ K” and 


The oxygen uptalce data obtained with the Uarburg 

respirometer, as described in section 4 ,2 ,3 • 1, is utilised to 
evaluate K':, K” and by Pujmoto^s Hethod. 


and iby i® plotted against y^, for the value of h as -> 

minutes. The ultimate oxygen uptalce of th& glucose solution • 
:&■£'! 'the 'dp4'b^n^ ,-turna’’'out to:'be 88G mg/1,-^ and ^ 

-it' "7 r ‘ '' .■ ' ‘ ?v'. • ' . t :;c 
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similarly from Pigux'S 'Zhf tho ultimate Oj-:yg€;n uptaJi.© 
of a micxX'bial suopensioii of Goncentration 300 mg/1, cornea out 
to be 395 mg/1, Therofoi'O/ the ultimate oxygen consumption of 

nij» Da./ wfl, 

the miorobio.l mass K" or dC/dx becornes! 1,3 2 tuii/w§^ The slope 
of this line gives Kq as 5.07 day) wliich corresponds to the 
autooxidation of the biomass at the rate of 93.7 ijercent per 
day/ vrhich is too high. Therefore for the purposes of the 
calculation of oxygen deficit by the rational model a value 
of Kq = 0.0514 day~^ and 0,1G2 day~^/ v/hich corresponds to 
the autooxldation of 5% and 15/* respictively biomass per day/ 
has been used. This value is more realistic as it has been ' 
reported by Sckenfelder (31). 


5,1.3 REAERATIOH R7I.TE C0K5T?4:iT Ko 

w ■ ’■ > " » >i.i. 

The data of the 0,0, concentration with time/ obtained 
as described in section 4, 2. 3. 3 are plotted as log (C< 3 ,-C) 
against time as a straight time in Figure 3 . The slope of 

this lines gives the reaeration rate constant K 2 as: 

Speed of the stirrer, rpm 30 ^0 70 

Reaeration Rate constant K 2 / day’”^, (base e) 4,146 7 .909 12.16 

5,2 SAG CURVES 

As describea in the section 4.1,4, the data for the 
D.O. profiles were obtained and plotted as b.O. deficit against 
time In Figures g to hi- for various initial microbial con- 
centration to substrate concentration ratio .The sag curves 
" as obtained from the Streeter and the rational 

rtipdel are also plotted in the same 


poted, that. for ail cases the yalue , of the^ critical 
by the ' Stre "■ Phelps first, order model is ■ 
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alv.'ayo lesn tViau th 
critical doficit is 
value. Tills will te 


e actually observed value nnd tlicit 
always greater Ui-an the actually 
ncl to predict higher critical D.O, 


of. the 
obser^'Gd 
deficit 


occuring earlier than tlio actual case. 


Hence any design -of a 


v/astev/ater treatment plant is likely to be over safe \vith 
tills model, as it is likely to suggest a higher degree of 
the treatment, than is required from the point of viev/ of 
the assivniliuting capacity of the strcarn for the pollution. 
Thus any design based on the Streeter - Phelps first order 
model is uneconomic, and a better appreodmation of the actual 
conditions can be acheived with the second order model propo- 


sed here. 

The sag curves plotted in Figure 2’ ■ to I'i- also shov/ 
that as the substrate to microorganisms ratio ie decreased 
v;e get more and more increase in D.O, deficit. This is obvi- 
ous because with an increased number of microorganisms work- 
ing on a substrate, the substrate v/ill naturally be metabol- 
ised faster and in process 3,0 .D. v;ill bo exerted at a higher 
rate, causing increased deplection of D»0, content* This dop 
endancy of the sag curve on .microbial concent rat icii also emp- 
hasises the importance of the microbial concentration for 
prediction of P,0, content in a water body . ^ 

The endogenous death rate as determined in the labo- 
ratory gave a value of K 0 higher than v/hat is usually encou- 
ntered in biological systems, as described above; in section 
5 ,>1. 2 . The values of the D ,6 * deficit from the rational model 
were caloAlated for values of to the^ 5 and 
15/4 antooxidation of th© micfobial mass per day, and are plO' 
tt ed in Figures g to 11 along with the , sag curve as Obtained 
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fro'in the Streeter - Phelpc modal and alfio the experimental data. 
Xt is seen* that the higher rate of antooxidaticvi gives a sag 
sarve nearer to the experimental valuec,, This means that in 
the stre.-iui models simulated in the laboratory the rate of 
endogenous respiration v/as higher tlian. tliat generally present 
in a ntreoiT\, On solving the equation 3.1.5 for higher voIligs 
oi a better fitting sag curve to the e:q 2 erin'iental data may 


be obtained. 

Similarly acijustmont in the values of ' all other const- 
ants v/ill produce a set a curves, vilaich may be tried and fitted 
to the observed sag curve in any stream and than tlio actual 
values of tVicae constants existing in the streojo may be appro- 
ximated wi'thout setting up an exx:>orimeatally simulated model 
of the stream in the laboratory. These constants then can be 
utilised in the design of a v;aste vjater 'treatment feiCility, 
discharging its effluent in the stre-im or can be put to 
some such other use. This gives us an indirect method of 
determination of Uiese constants £oz' a stream. This also 
adds to the flexibility of the model as there are seven of | 

them and variation in every one will change the sag curve. 
Therefore, by trial and erxx>r almost any data observed in } 

the nature can be filled to a curve and the constants deter- ; 

mined for the prediction of D.O. downstreaiTi to v/astev/ater i 

out fall ♦ ' ; . ■ P . ■ , . ; , ■ -r ■ -fi 

In figures and 9 v/e see that the p «0 • deficit . | 

decreases a :littl^/ in the beginning, as predicited by the | 

■^rational' model . ' This h® ^dae5::to ':tbe;'higher-^rate' y::| 
erationy due to the high yal'ue of : d beginning / . ^ 




vjl'iich it: iaoj.o than bhe rat© oi dooiiy^enation thus offsetting 
it. A© tiro© passes the D,0« deficit decreases and i/ith 
it., the rat© of ruaer-atiou decreases too. A tiw© comos wlieti 
tbo iotCi of reaGration decreases to the ei'tsnt so as to 
become lov/er than the rate of deoxygonation, tad the D.O. 
deficit again starts increasing. The D.O. deficit starts 
roduciug once again, vhen . the rate of deoxygensitrion rodaces/ 
due to reduction in rosniainiirg B.O.D., sufficiently, so as 
to become lovior thcan the rate of reaeration. Th© ©Kpeilrnental 
data docs not shovj the initial decreases in D.O. deficit. 

This could be due to a possibility of a mistake in measurement 
of tl'fe volume of stock suspension to produce such lov/ rnicro~ 
bial concentrations* 

The sag cura/es as plotted in figure 3^ to Figure 1 2- show 
that for smaller values of t, the values of th© oxygen deficit 


as obtained from the Streeter - Phelps Model and those from 
the rational Model vary significantly. But as is clear .from 
Figure 12- , for higher values of initial microorganisms 

concentration, (Xq~ 10*5 mg/1) in the presence of lov/er sub- 
strate concentration (higher values of t) the tv/o profiles 
agree with each other to a significant extent. This is in 
contrast to the fact that the Streeter - Phelps model dis- 
regards the microbial concentration while predicting the D .0 . 
Actually as the tirne elapses, more and more substrate gets 
converted to microb-ial mass and therefore, for higher values 
of t, the time, the effect of the miGrobial cohcent ration 
should be more pronounged . Therefore it can reasonabiy be 

data frop the two models should differ /i 
significantiy for smaller values of t and microbial congeh'*' 

the .time, passes,- - substrate: concent ration 



CiCuJ- UUtA LilCj llliCLOijXClX COUO'iliCfQVtiOli Ut:iCX'<Sc»iJti!3 • i'J 

XllviC C.UO JjUi.)i;bX'uL»i: COliC’^lVLXiUt ioll btiCQUitstJ SO 
LiiliiL iL cuwou uo luiva a uiyuiXicatic tiXlocX/ ciuci uiit- 
ci^i^ouLxaiXy ijocoiii^iLi j. iiixs'c ordt^r ouw# dtiptiudiug ouly 
Liie micxobiai ouucentxatiou. Tims ultiiough the Streeter - 
i-’iieips toriiiuiatiou is suhsti'ate concent raciou dependent 
and tlie ratioudi iTodei becomes microbial concentration 
dopoudeut rox the higher values ol t, both the expression*^ 
are tlie IxrsL order ones and hence shov*/ sifnilar values* 

A;, lar as Uio actual values oi. tlie oxygen deliti^ 
are concerned we note that the experiiviental values are 

liigher tiiau those predicted by tiio rational model and are 

lower in the bcyinuiuy and hiyher later than those yi'^'ti* 
by tiie btreoter r-licips Model , Tliis Ciperimental data 
also closer to tiie rational Model chan the Streeter 

• fc 1 Y\ 

lormulation, ns' a vaatcer ol fact the actual data rs 
ayrecnieut with the btreciicr - fc'lielxjs valus only in a ve "V 
small iusiyuiricant region v/here the two plots cross ea 

oUier* ht all other times the experimental values tally 

» ■ m *1 ■ ■ 

better with tlie rational Model rather than the clas^>xt;' 

~^■r •^ry 

Streeter - Phelps first order formulation/ which a/. 

■■ 

^ ^ ^ she results predicted by the Streeter 

■ ... T ■ 

Model the critical couctition^ jxb predicted by c it. 
are distributed over a relatively long «duration o± 

weii 


a sort of flat region irl the curve* Tl^his again 9 

4 . to the impo rt afice b i« the niicrohial concent ration * 
period of time enoygh biomass is produced to hayt- n 

rate of respiration of its own* fvt tlie same tlm 
■bolder and dead rnicroorgahisms to be 


•igni 




fict 


iciucing 

se;;.' y; 

ain 
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ooii v ^ ioti Ciits 1 X xlico ilii o X'O i^i iXX itim.;/i^ cj.o^l» 

iiO C XC^ OXXciOC Oti CiX*^ X^ *0 * OOiiLc^H'ti OX CiiO Ot^XOcUii/ 

hath axti yoii'iy to cause deoxygenatioii and lionctx uiaiiitaxa its 
x'uto , ii uaxance xs xeacliou. oetwuen ciit xute ox cleojiygeuatiou 
and x'eoixygeuation and the D«0, remaius coiistaiit at that level 
xo X ^onio cxiiiO « niitii enougii suoscxaee aixo iuxcxox)xal ttiass has 
been completely o>:idised to G 02 and vyatex/ the phase of 
recovery strts and the rate of deo;xyyeuation falls belovy 
che x’ate or reoiiyyeu.iition causing the b.O. deficit to rise ayaiii 

The dtreater - t^helps equation, as it does liot give aviy 
covisideration to tiie coviceiitration of the micx’otiai mass, gives 
tlie scuae sag curve for all d/X ratios* It is clearly shown in 
Figures 8 to 12 that tliis is not so. The Increasing vuicrobial 
concentration causes higher D.O* deficits occuriug earlier 
than in the case of lov/er microbial concentracious* This is 
also cox'roborated by the sag curves predicted by the ration- 
al luodel . • 

ns the constants were determined at a different conce- 
ntration of substrate end the microbial mass, it is quite likely 
that their values for the simulated streara arc not the same as 
in the set up used for their detexmiiiation. This perhaps 
accounts for the poor correlation of the ejcperiraental data 
with the rational model sag curves* 





COi'jCiiU iJXOii> u 


foliowiay conclusion a are urawn iroui •Lhia a'cudy; 

1. '1‘he a’creeter Phelps iuoilcl in noc adequate for the 
preciiction of the sa 9 curv'es dua to the foiiowinq raasono: 

(a) The, critical conditions at predicted by this itiodel 
occur earlier and are jiiore severe, than those obser'ved experimeu- 
tally, 'iiiic v;ill tend to produce oversafe designs for wastewater 
treatvacnt plants. and hence will not be econoniicaliy desirable. 

(b) This iiiodei predicts that the sag curve revnaine the 

saiue for the same intial substrate concentration irrespective of 
the microbial concentration present in a strearii. dxpex'iriicntal data 
shows that it is not so. As the iidcrobial covicentrations along the 
couse of a stream are apt to vary significantly, depending upon 
the. runoff and wastewater discharge received by it, any prediction 
of the b.O. deficit by the Streeter - Phelps model are. likely to 
be erraneous. This discrepancy is taken care of in the rational 
model as different sag curves are predicted for different initial i 
microbial concentrations. , j 

2, The constants determined under conditions different than ' 
those in the simulated stream liave caused the poor correlation of the; 
oicperfmental data to the sag curve predicted by the rational model. | 
However tVie rational raodei: still gives a be-te.ter approximation to r 


the experinientai data, as corapax'ed to Streeter Phelps rnodel, show- 
ing that the inetabolic process of the pollutants is a. second 
order are rather than of first -order. By a judieions choice of " 


CGnstaiits it will be possible to acheive a significant cofr elation; 
in the experimental and' rational model sag :ai ryes . ^ > 
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i\PPBHDIX .. 



*A‘ 


1 IllTEGRATlOM OF SQUATIon 3,6 


Equation 3,6 is 


raa:c 


6CA~3) 
Kr,+ S 


Sepernting the variables wg get/ 


5 

“5(a..S) 


dS 


■" Hma:*: 


• r dt 


* « , A , X • 1 


oi/ A 


• + (1+ K-j^/A) , 


1 


] 


dS = 




integrating from S == 3^^ at t = 0 to S - S at t -- t vjo get 


^ In (S/So) - (It K^/A) In 


^ * * A«l»2 


Therefore, t 


1 

l^max 


*■ I 'Si,'*"' ! 


or. 


t 


u, 


In 


l^A~So * S 


max 




.'A 


^ y 

"^o j 


. A*X.3 
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APPKMDIX ‘ B ' 


B . 1 i':vAl.LJATji:ow o?-‘ f::oi:TST.M'^TS 

B*l»l Kq, 


Time, t Microbisl concentration, x Sutntrate 

hours ^'^3/3- tiori as ^ 

nvg/1 


0 

Li 

440 

5 

12 

430 

9 

13.5 

430 

13 

13 » 3 

430 

16 

20 

400 

24 

27.5 

370 

27 

38 

350 

30 

55 

2.10 

32.5 

63 

330 

36 

84 

250 

39 

94 

180 

42 

107 

140 

48 

115 

130 

51 

122 

no 

54 

134 

100 

57 

140 

90 

60 

150 

100 


153 

30 

12 

.158 ■ 

i/ . 30 


concent ro' 
0 « ij ♦ 



B, 1.2 DETaRflliNJAriOIl OP K*' AMB V.r^ OM VmR'iUr Yit 
liSSPIHOMBTSK 


B-.1..2.1 D3Ti!;i'!MIN/VriOB 0:F ELABK COM?3'rA^:TS 


Teroperaturo 


■D 

*• o 


V 


f 


5= 25'^U 

- lO',QOO mm o.t Bradie's fluid 
™ 0.0233 

=5 2.6 ml - 2600 fil 


Flash Mo . 

Vg, ml 

FlasJc constant K , l/rnrn 

1 

14.2 

1 • o6 

2 

13 «u6 

2.06 

3 

16.33 

1,79 

Kinco, 

X 

K Kh in jiAl of ojcygon consamod nt 

MTP 

32 


■— Kh in mg/1 of substrate 
22,4V 


wliore. 


0.53 Kh 


=! volumn Gubstrate, ml 


V 



♦ 


63 

B. 1,2.2 

OXYOBN UFTJ\KI2 DATA 


Time, t 
minutes 

Oxygen uptalc© for glucose 
mg/1 

O^'ygen uptotinlie for 
mi croo rga n i sm anuo r 
c n do g e nous co d di t i o n , 
rng/1 

30 

46.3 

12.6 

60 

103.0 

36,5 

90 

166*2 

73*0 

120 

220.0 

100.0 

160 

254.0 

115,0 

IBO 

272.0 

124,5 

210 

301.0 

139.0 

240 

322.0 

151.0 

270 

342 .0 

151.0 

300 

366.0 

176.0 

330 

334.0 

188,0 

360 

405 .0 

214.0 

390 

432.0 

228.5 

420 

455.0 

246.0 

450 

476.0 

262,5 

430 

494.0 

273.5 

510 

506,0 

284.0 

540 

542.0 

300.0 



s.1.3 Mh 

LUATIOH OF Ko , Till 

RF«SIidTIOM RATO 

* 0 ^.^ f'- 


T ernxj e r a tu r o 

- 22 ^'^C 




0,8 mg /1 



Vo I , of v/a tor 

= 2.5 1 



Vessel used 

; 30 ern, :: lO era 
D , 0 , concent rat 

i4 cylindrical Ja. 
.ion, rag/l 

TIMS, jrdnut 

a X . 30. r pm X 

bo rprn. X 

70 rprr. 

0 

1 • '9 4 

r* 9 "7 

5 .00 

10 

2.18 

3 .54 

5 , 47 

20 

2.35 

0,80 

5,63 

30 

2.57 

. 10 

5.67 

40 

2.89 

4.25 

6.10 

50 

3.05 

4- ,50 

6.27 

GO 

_ 

4.02 

•5,45 

70 


5.00 

6 ,70 




; O' 


B « 2 DJVTA Foa gjRVES 

B. 2.1 E;TI2U1MBKTAL VAIiUBS 


Time, hrc. X 


Oxygen Deficit, mg/1 


XSot No.l Xset Tlo.2 Xsot M0.3X Get no.4i) Set !'o. 5 


X 3:q=2 , 1 


X 2 ;' s "4 ,2 X . 3 X Xq= 3 . 4 X :ro-lO . G 


^ mcr/1 X_rria/1 X rn'j/l X rng/l X rnq/I- 


0 


0.9 


1 

O 

• 

o 

2 

1.0 

3 

1.0 

4 

1,0 

5 

1.10 

6 • 

1.10 

7 

1.0 

8 

1 , 10 

9 

1.2 

10 

1.2 

12 

1.2 

16 

1.2 

20 

1.2 

24 

0 . 

o 

28 

1.10 

32 

0.90 

3G 

0 .90 

40 

0 , 8 

44 

0.8 

48 

0.6 


0.9 

0.9 

0,95 

0,95 

0.95 

0.95 

1.0 

1.0 • 

1,1 

1.0 

1.1 

1.1 

1.1 

1.15 

1.15 

1.20 

1.2 

1.2 

1.2 

1.3 

1.25 

1.3 

1.2 

1.25 

1.25 

1.25 

1.3 

1.2 

1.3 

1.1 

1.2 

1.1 

1.1 

O • 95 

1 «o 

0,90 

1.0 

0 ».8 

0.9 

: 0..8 , 

0.8 

0.8 


0,9 

0,9 

1 .0 

1.0 

1.0 

1.2 

1.1 

1.4 

1.15 

1.4 

1.15 

1.5 

1.1 

1.6 

1.2 

i.Gi; 

1.25 

1.7 

1.3 

1.7 

1.3 

1.7 

1.4 

1.6 

1.4 

1.6 

1-.3 

1 .5 

1.2 

1.5 

1.15 

1.3 

1.0 


1.0 

1.1 

1.0 

.0,9 

0.9 

0,8 

0,9 

0.7 


3.2.2 
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SAG CURVE BY ^.^^TXOfs!AL MODEL 


a, 


2.1 V/\LUSS rOR ~ 0.051‘i 


X 

3!:„i=2.1 mg/1 X 

4,2 rng/l X 

7} 

u * J Ti 

g/1 ’X 

3.4 m 

g/1 ^ X 

lO.'sTn 

^g/1 

3 I 

hi.-.'S .X 

a X 

mg/1 1 

t I D~l'~ 

ht's Irag/Xl 

t X D X 
hrs Xmg/1 1 

■'T^T 

hrr: X 

D X 
rng/l X 

" I 

hrs X 

D 

mg/l 

15*0 

0 

0.9 

0 0.9 

o 

0.9 

o 

0.9 

0 

0.9 ■ 

:4.o 

4.57 

0.456 

2,33 0.324 

1.61 

1 .012 

1,21 

1.1651 

0 *9^7 b 

1.2G5 

:3 ^0 

3*65 

0.4-52 

4*6G 0*320 

3 .19 

1 .023 

2.42 

1*2 uuJ 

1.96 

1 ♦ 4 ci 

:2.o 

12.37 

0.480 

0 # Q 6 0 • B 3 1 

4.75 

X ♦O.j o 

3.54 

1*2 cio3 

2,95 

1.5325 

:i.o 

15 .63 

0*511 

9.00 0.333 

5,31 

1 .0441 

4 « 85 

1,2301 

3 ,95 

1.5317 

O 

* 

p 

19.09 

0 « 539 

11*10 0*3^b 

7.33 

1 *0442 

6 *08 

1 , 287 

a- e 9 C 

1 . 5307 

.9.0 

21*19 

0 . 562 

13*17 0*836 

9 .40 

1 .0440 

7 .32 

1.286 

5,99. 

1.5293 

.3.0 

25.19 

0.500 

15.22 0.337 

10.96 

1.0435 

3.5B 

1.285 

7 ,05 

1.527 

.7.0 

23.10 

0.592 

17.26 0.337 

12.54 

1 *0426 

2.86 

1.284 

8.31 

1.525 

l5*0 

30.97 

0,6o4- 

19.32 0.83G 

14.14 

1.0414 

11.17 

1.262 

9.24 

1.523 

l5,0 

33.81 

0,610 

21.39 0.835 

15.77 

1 .039 

12.52 

1.279 

10.39 

1.520 

L4,0 

36.65 

0.6099 

23.51 0.334 

17.45 

1.037 

13.91 

1.276 

11,58 

1.517 : - i 

L3.0 

39.52 

0.609 

25.67 0.332 

19 . 19 

1.034 

15.39 

1.273 

12,62 

1,513 1 

L2.0 

42.43 

0.600 

27.90 0.329 

2.1.00 

1 .0<i> 1 

16. SB 

1.2G9 

14.13 

1.500 " 1 

LI .0 

45.43 

0.605 

30,22 0.826 

22.89 

1.026 

13.48 

1.264 

15.52 

1.502 . I 

10 ,0 

^i‘8 . 5^ 

0.602 

32.66 0.306 

24.89 

0.992 

20.18 

1.138 

17 ,00 

1.399 - ; 1 

1 

9.0 

51,32 

0 « 592 

35,25 0.757 

27,04 

6,923 

22.01 

1 . 105 

13.59 

1.289 j 

8.0 

55.29 

0.556 

38.03 0.704 

29.35 

0,850 

23 .99 

1.0G7 

20.32 

1.165 .. i. 

7.0 

59 .06 

0.511 

41,06 0.640 

31,39 

0.773 

26,18 

6,909 

22.24 

1.051 : j 

6.0 

63.21 

0 . 464 

44,44 0,576 

34.74 

0,693 

28.6 

0,813 

24.4 

0,9373 : [: 

5.0 

67.92 

0.412 

48,29 0.509 

38,0 

0.692 

31,47 

0.712 

26,90 

o.sm -i 

4.0 

73.45 

0.355 

52,85 0.435 

31.88 

0.518 

34,84- 

0,604 

29.89 

,;0 ,692 ; ■ ; ,jj 

3.0 

80.29 

0.291 

58,58 0.355 

36.73 

0.420 

39,08 

0.488 

33,65 

■r 0 . 557 

'■ 

• • ‘ }'■ 

2 .0 

89.57 

0.220 

66.28 0.266 

53.38 

0.314 

44.90 

0,362 

: 38.S3:- 

; 0 ,* 



t,2.2 V/iLUKS ?0H K. 


0.162 


SUvOctrat 
concent r 
tion. s. 

aX x-i - ; 

2.1 mcj/li; 

^^0=^4* 2 mg/^1 1 .3 

mg/1 X:v 3 :^ 8.4 

rng/1 X 

b-10.5 

' m/i 

a* t j 

1 hrs ) 

i mg/1 1 

t 1 

lirs 0 

13 X t } D X 

rng/1 X lire X rng/1 X 

t X 
Xrrn X 

o nr 

^g/3- X 

lijr's X 

yj 

rng/'X 

rng/1 

J 










2 5.0 

0 

0.9 

0 

0.9 

Q 

0.9 

0 

0.9 

0 

0.9 

2 4 .0 

4.57 

0.476 

2.38 

0.859 

1.61 

1,055 

1.22 

1.20 

0,98 

1.305 

23,0 

8 .65 

0.475 

4 .66 

0 a 86 

3,19 

1 .087 

2.42 

1.371 

1.96 

1*552 

2 2 .0 

12.37 

0 .508 

C ,86 

0,867 

4.75 

1 .099 

0 #o4 

1,371 

2.95 

1 . 620 

21.0 

15 .83 

0 . 546 

9.00 

0.369 

6.31 

1,110 

4' » 3 L> 

1*8 i 1 

3 .9 5 

1.619: 

20.0 

19.09 


11.10 

0,8/1 

7.65 

1.121 

G.08 

1.370 

4,96 

1.619i 

19 ,0 

22.10 

0 . 600 

13 . 17 

0.872 

9.40 

1.121 

i ^ >lj AL 

1.370 

5.99 

1.6l7j 

IB.O 

25,19 

O.C24 

15.22 

0.874 

10.96 

1.121 

8 • 58 

1.369 

7,05 

i.eiei 

17.0 

28,10 

0.643 

VI *26 

0.B82 

12.54 

1,121 

9.86 

1.368 

B.13 

l*6i4j 

16.0 

30.9 

0.656 

19.3 2 

0 ,36 1 

14.14 

1.119 

11.17 

1.366 

9.24 

1.612 

15.0 

33.81 

0.663 

21.39 

0 .8o 1 

15.77 

1.117 

12.52 

1.364 

10.39 

1.6 id 

i 

i 

14 ,0 

36.65 

0.664 

23.51 

0.880 

17 .45 

1.115 

13.91 

1.362 

11.58 

1.60l| 

13.0 

39.52 

0 . 663 

25.67 

0.878 

19.19 

1.113 

15.36 

1.358 

12.82 

1 . 6o:| 

i 

12.0 

42.43 

0.662 

27 .90 

0.876 

21.00 

1.110 

16.88 

1.355 

14.13 

1.599 

11.0 

45 . 43 

0.661 

30,22 

0,874 

22.89 

1.106 

18 « 48 

1.350 

15.52 

■ f 

1.59f 

10.0 

48.45 

0 , 658 

32.66 

0.369 

24 .89 

1.101 

20.18 

1.344 

17.00 

! 

■1.57^ 

9.0 

51.81 

0.654 

35,25 

0.864 

27.04 

1.054 

22.01 

1.252 

18.59 

1.46^ 

8.0 

55.29 

0,643 

38 <.03 

0,811 

29.35 

0 ,989 

23.99 

1.164 

20,32 

1.34| 

7.0 

59.06 

0.604 

41.06 

0.753 

31,89 

0.910 

26.13 

1.069 

22.24 

i.22j; 

■ ■ 1 

6.0 

63.2 

0.559 

44,44 . 

0.692 

34.74 

0 .826 

28,64 

0.964 

24.40 

1 i id 

1 

'5,0 

67.92 

0.506 

48,29 

0.621 

38 ,00 

0.741 

31.47 

0.864 

2,6.90 

o.m 

:■■■ I- 

4.0 

73.45 

0.451 

52,85 

0.550 

41.88 

0,653 

34.84 : 

0.758 

29 .89 

:■ / Iv 

3,0 

80.29 

0.391 

53.53 

0.474 

46.73 

0.559 

39.08 

''.0,646 

33.65 

■0.*7l 

CO 

O 

89,57 

0.323 

66,28 

0,389 

53 ,38 

0.456 

•44*90 

0,525 

:38’.83; 

^rb/sl 






B, 2«3. SAG 


BY STRJiBTSI'i ~ PHBhPS T-iODBL 



= 0.04606 hr"*^ 



= 12.16 day’”^ « 0.507 hr*"^ 



= 0.9 rng/1 

Time, t 
hours 


D.O. Deficit 
^ 9/1 

0 


0.9 



1.65b 

4,03 


1,77 b (oritical values) 

B ' 


1.496 

12 


1.262 

16 


1,050 

20 


0,370 

24 


0,729 , 

23 


0 * 6o5 • 

32 


O.bOb 

36 


0.420 

40 


0.349 

44 


0.'292 

. 48 


0.241 


APPENDIX 'C' 


fORlRAM I V ( IBm 704A/1401 VERSI0N> PROGRAM FOR TriE SOLUTION OF 1 HE 
Equation 3 . i .5 » sy runge-kutta ( fourth oruer ) method. 


•DIMENSION AKAY< 10) 

, ,, F.'(;S >,D.,)>'C * ( AKS'TS ),*D/ ( S^F( A-S ) > - I AkDAS+B* { AKS+5 ) /P 1 
./T|:;ffl|i^EAOy:T.»'AKX >AkS »AKMAX:,AK 2 , AiCDAS,AKt ,DCDx 

ao;i .f:;QAtT;7Fio!.4) : . 

Q u VPl||/S|if|-3T^i<XTAK S , AKMAX , AK2 , A KDAS , AKE , BCDX 

.3»5X»3Hk 5 = E8.3 > 5X » 5HKMAX = f q . 4 /X >3n:<2 = Ea.3>3X» 5 HKp 
V . I V7 X , 3 H X E = F 8 . 3 , 5 X , 6 H D C / D X = F 8 . 3 ) 

V ■ ;;ul i'o 2 ' » D E L » E P S , I T E R » J S E T S 

• ibZ F0RMAT(2F8,3,214) 

PRINTIIA, DEL, EPS, ITER, JSETS „ - ,..,,Errs r;v 

ITA Format (X, 4 HDE|. = F8.3»5X.4HEPS = F8.3 »5X,iqHITtiRATl0NS AlLOWEu= I 4 , 5X , 5 
iHSETS=lA): 

6= ( AkE4:-DCDX ) / ( AkX^AKMAX ) 

C=AR2/AKMAX 
DObL=l»J3tT5 
PRINT 103 
printioa 
. PP I NT 103 

( X > 1 4 ( 1H“ ) ) 

IAHNEW SET BEGINS) 
v ^,;X0. 5 DO > H » s M I N 
:-'l''D5''Ed(^M'AT(5;F0V3) 

A-AKX^fXO + SO 

PRlNT106»SO»XO»DO,H>SMlN ^ ^ ov rLJc•^.fwi 

1U6 FnRfAT(X»3HS0 = F8.3»3X,3HX0 = F8»3»3X , 3HD0 = F 8 *3 » 3X , 2HFtsF8 . 3 » 3X » SFtSMlN.s 
i==F8.3) " 

INTERS { SO-SMIN )/( -DED + l , ■■ - .:r -T.-. ■•vQ!' 

D=D0 - ■ ■ ■ - • ' 

PRiNTloT » INTER ... ' . 'Q- yU-?::V®S 

1C7 FOFMAT(X»10HINTERVALS=I6)' ■ r',,’. ■ ' ’ ' • , 'Fp-TV 

i' PRINT112 ; ■ V 











KilSIi 


fvjlK*'! •;?0 •' •;T|Fr -F, 
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iJ9 i"ORMAF (3X»1HI ,4X ,9HS IN MG/L » 4X » 1 3HT 1 Mt I fv HOUR 3 » 1 X > 1 5 riOEF I C I T IN 
1MG/L,1X»16HERR0R IN DEFICIT) 

112 F0RMAT(X»65'{1H*)/ ) 

DO 7 I*;l, INTER . 

» D ) 

S+H/2 . , D+A 

^;'pF£;V^i<«T4FH*F(S+H»D+AXAY3 ) 

C0RR=(AXAY1+2.'«-( AKAY2+At<;AY3)+AKAY4) /6, 

DNt.W”D+CORK 
IF (DREW) i0,9,9 
10 PR1NT113 

115 FOFU-iAT(X»311-IOXYGEN deficit BECOMES NEGATlVt) 

9 ERRORS I DNEW-D ) /DREW 

IF(AB5(ERROR>-EPS)ll»li,3 

3 r --DNCV/ 

5 .CJNTINUE 

(] A L C U L A T I 0 N 3 F OR ' T ' » G I V t N * S * 


11 AS=A-S 
• ■ AS0=A-50 

cc'V -..xi-= alog ( so ) “ALO g ( s ) 

{ as ) -alog i aso ) 

• ) -iJ-XZ-v- ( /M<3/ A ) irXl ) / /vKMAX ) 
T » DNFW » ERROR 

E 1 3 . 6 , 2 X ) ) 

:\iv:''''7;’ CONTINUE 
■ A 8' C(j#iNUE . 

STOP 
END ■ 




(V£ - - 
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